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ABSTRACT: In this paper, we introduce an approach to increase density of field-effect transistors 
framework instrumentation amplifier input bias circuitry. Framework the approach, we consider the 
manufacturing of inverter in heterostructure with specific configuration. Several required areas of the 
heterostructure should be doped by diffusion or ion implantation. Moreover, do pant and radiation defects 
should by annealed framework by optimized scheme. We also consider an approach to decrease value of 
mismatch-induced stress in the considered heterostructure. We introduce an analytical approach to 
analyze mass and heat transport in heterostructures during manufacturing of integrated circuits with 
account mismatch-induced stress. 
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INTRODUCTION 

At present, several problems of the solid state 
electronics (increasing performance, reliability, 
density of elements of integrated circuits: diodes, 
field-effect and bipolar transistors) are intensively 
solving [1-6]. To increase the performance of these 
devices it is interesting to determine the materials 
with higher values of charge carrier mobility [7- 
10]. One way to decrease the dimensions of 
elements of integrated circuits is to manufacturing 
them in thin film heterostructures [3-5,11]. In this 
case, it is possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12] and 
development of epitaxial technology to improve 
these materials (including analysis of mismatch 
induced stress) [13-15]. An alternative approach 
to increase the dimensions of integrated circuits 
are using of laser and microwave types of 
annealing [16-18]. 


In this article, we introduce an approach to 
manufacture the field-effect transistors. This 
approach gives a possibility to decrease their 
dimensions with increasing their density 
framework, an instrumentation amplifier input 
bias circuitry. We also consider the possibility to 
decrease the mismatch-induced stress to decrease 
the quantity of defects, generated due to the stress. 
We also consider the heterostructure, which 
consist of a substrate and an epitaxial layer (see 
Fig. 1), a buffer layer between the substrate and 
the epitaxial layer. The epitaxial layer includes 
several sections, which were manufactured by 
using other materials. These sections have been 
doped by diffusion or ion implantation to 
manufacture the required types of conductivity (p 
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or n). These areas became sources, drains and 
gates (Fig.1). After doping it is required to 
annealing of dopant and/or radiation defects. The 
aim of this paper is to analysis of redistribution of 
dopant and radiation defects to determine 
conditions, which correspond to decreasing of 
elements of the considered filter and at the same 
time to increase their density. At the same time, 
we consider a possibility to decrease mismatch- 
induced stress. 


Fig. 1b. Heterostructure with a substrate, epitaxial 
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METHOD OF SOLUTION 


To solve our aim we determine and analyzed spatio-temporal distribution of concentration of dopant in the 
considered heterostructure. We also determine the distribution by solving the second Fick's law in the 
following form [1,19-23] 


6C(x,y,zt)_ 6 pola y,2t) 4 poe lay.zt) 9 po luy.zt) 


Ot Ox Ox Oy Oy OZ OZ 
0| D, ty 
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Q—|—V »y,Z,t)| C(x, y,W,t)dw 
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with boundary and initial conditions 
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Here C(x,y,z,t) is the spatio-temporal distribution of concentration of dopant; © is the atomic volume of 
L 

dopant; Vs is the symbol of surficial gradient; [Cc (x, y, Z, t)d Z is the surficial concentration of dopant on 
0 


interface between layers of heterostructure (in this situation we assume, that Z-axis is perpendicular to 
interface between layers of heterostructure); /“(x,y,z,t) is the chemical potential due to the presence of 
mismatch-induced stress; D and Ds are the coefficients of volumetric and surficial diffusions. The values of 
dopant diffusions coefficients depends on properties of materials of heterostructure, speed of heating and 
cooling of materials during annealing and spatio-temporal distribution of concentration of dopant. 
Dependences of dopant diffusions coefficients on parameters could be approximated by the following 
relations [21-23] 


V"(x, y, zt) 


vy | 


af SY, Zot V(x, y, Zt 
D, =Dy,(x, y,z.T)) 1+ eee 16 HO ge .(2) 


Here Dz (x,y,z,T) and D1s (x,y,z,T) are the spatial (due to accounting all layers of heterostruicture) and 
temperature (due to Arrhenius law) dependences of dopant diffusion coefficients; T is the temperature of 
annealing; P (x,y,z,T) is the limit of solubility of dopant; parameter y depends on properties of materials and 
could be integer in the following interval y €[1,3] [21]; V (xy,z,t) is the spatio-temporal distribution of 
concentration of radiation vacancies; V* is the equilibrium distribution of vacancies. The concentrational 
dependence of dopant diffusion coefficient has been described in details in [21]. The spatio-temporal 
distributions of concentration of point radiation defects have been determined by solving the following 
system of equations [19,22,23] 
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I (xy,z,0)=fi (x),2), V (xy,2,0)=fv (%y,2). 
Here / (x,y,z,t) is the spatio-temporal distribution of concentration of radiation interstitials; I’ is the 
equilibrium distribution of interstitials; Dj(x,y,z,T), Dv(x,y,z,T), Dis(%y, Z,T), Dvs(x,y,Z,T) are the coefficients of 
volumetric and surficial diffusions of interstitials and vacancies, respectively; terms V2(x,y,z,t) and J2(x,y,z,t) 
correspond to generation of divacancies and diinterstitials, respectively (see, for example, [23] and 
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appropriate references in this book); kiv(x,y,z,T), kii(xy,z,T) and kvv(xy,z,T) are the parameters of 
recombination of point radiation defects and generation of their complexes. 


Spatio-temporal distributions of divacancies @y (x,y,z,t) and diinterstitials ® (x,y,z, t) could be determined 
by solving the following system of equations [19,22,23] 
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D (x,y,2,0)=far (x,y,2Z), Pv (x,y,2,0)=fov (x,y,Z). 
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Here Dai(x,y,z,T), Dov(x,y,z,T), Dais (x,y,Z,T) and Dovs(x,y,z,T) are the coefficients of volumetric and surficial 
diffusions of complexes of radiation defects; ki(x,y,z,T) and kv(x,y,z,T) are the parameters of decay of 
complexes of radiation defects. 


Chemical potential sm in Eq.(1) could be determine by the following relation [19] 
M=E(Z)Q oj [ui(X,y,2,0) + uji(x,y,Z,0)]/2, (7) 
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is the deformation 


tensor; ui, uj are the components u:(x,y,z,t), u(x,y,z,t) and u-(x,,z,t) of the displacement vector 
u (x, y,Z, t); Xi, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 
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where cis Poisson coefficient; & = (as-az1)/azz is the mismatch parameter; ds, azz are lattice distances of the 
substrate and the epitaxial layer; K is the modulus of uniform compression; / is the coefficient of thermal 
expansion; 7; is the equilibrium temperature, which coincide (for our case) with room temperature. 
Components of displacement vector could be obtained by solution of the following equations [25] 
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heterostructure, dj Is the Kronecker symbol. With account the relation for oj last system of equation could 
be written as 
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We determine spatio-temporal distributions of concentrations of dopant and radiati-on defects by solving 
the Eqs.(1), (3) and (5) framework standard method of averaging of function corrections [26]. Previously 
we transform the Eqs.(1), (3) and (5) to the following form with account initial distributions of the 
considered concentrations 
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Further we replace concentrations of dopant and radiation defects in right sides of Eqs. (1a), (3a) and (5a), 
are not yet known average values aip. In this situation we obtain equations for the first-order 
approximations of the required concentrations in the following form 
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Integration of the left and right sides of the Eqs. (1b), (3b) and (5b) on time gives us possibility to obtain 
relations for above approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and radiation 
defects by the following standard relation [26] 
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We determine approximations of the second and higher orders of concentrations of dopant and radiation 
defects framework standard iterative procedure of method of averaging of function corrections [26]. 
Framework this procedure to determine approximations of the n-th order of concentrations of dopant and 
radiation defects we replace the required concentrations in the Eqs. (1c), (3c), (5c) on the following sum 
Onpt P n-1 (X,y,z,t). The replacement leads to the following transformation of the appropriate equations 
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Integration on the left and the right sides of the Eqs. (1d), (3d) and (5d) gives us possibility to obtain 
relations for the required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the following 
standard relation [26] 
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Further we determine solutions of Eqs.(8), i.e. components of displacement vector. To determine the first- 
order approximations of the considered components framework method of averaging of function 
corrections we replace the required functions in the right sides of the equations by their not yet known 
average values ai. The substitution leads to the following result 
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Approximations of the second and higher orders of components of displacement vector could be 
determined by using standard replacement of the required components on the following sums aj+u((x,y,z,t) 
[26]. The replacement leads to the following result 
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Framework this paper, we _ determine 
concentration of dopant, concentrations of 
radiation defects and components of displacement 
vector by using the second-order approximation 
framework method of averaging of function 
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DISCUSSION 

In this section, we have analyzed the dynamics of 
redistributions of dopant and radiation defects 
during annealing and under influence of 
mismatch-induced stress. Typical distributions of 
concentrations of dopant in heterostructures are 
presented on Figs. 2 and 3 for diffusion and ion 
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corrections. This approximation is usually enough 
good approximation to make qualitative analysis 
and to obtain some quantitative results. All 
obtained results have been’ checked by 
comparison with results of numerical simulations. 


types of doping, respectively. These distributions 
have been calculated for the case, when value of 
dopant diffusion coefficient in the epitaxial layer is 
larger, than in the substrate. The figures show, that 
inhomogeneity of heterostructure gives us 
possibility to increase sharpness of p-n- junctions. 
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At the same time one can find increasing 
homogeneity of dopant distribution in doped part 
of epitaxial layer. Increasing of sharpness of p-n- 
junction gives us possibility to decrease switching 
time. The second effect leads to decreasing local 
heating of materials during functioning of p-n- 
junction or decreasing of dimensions of the p-n- 
junction for fixed maximal value of local overheat. 
However framework’ this approach of 
manufacturing of bipolar transistor it is necessary 
to optimize annealing of dopant and/or radiation 
defects. Reason of this optimization is following. If 
annealing time is small, the : did not achieve 
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any interfaces between materials of 
heterostructure. In this situation one cannot find 
any modifications of distribution of concentration 
of dopant. If annealing time is large, distribution of 
concentration of dopant is too homogenous. We 
optimize annealing time framework recently 
introduces approach [15,25-32]. Framework this 
criterion we approximate real distribution of 
concentration of dopant by step-wise function (see 
Figs. 4 and 5). Farther we determine optimal 
values of annealing time by minimization of the 
following mean-squared error 


a7 [file (x, y,z,O)-w(x, y,z)]d zd ydx, (15) 
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Fig.2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Increasing of number of curve corresponds to 
increasing of difference between values of dopant diffusion coefficient in layers of heterostructure under 
condition, when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant 


diffusion coefficient in substrate 
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Fig.3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, which 
is perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to annealing 
time © = 0.0048(Lx2+L,2+L22)/Do. Curves 2 and 4 corresponds to annealing time © = 0.0057(L.2+Ly2+Lz?) /Do. 
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urves 1 an corresponds to homogenous sample. Curves 5 an corresponds to heterostructure under 


condition, when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant 


diffusion coefficient in substrate 
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Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 


Increasing of number of curve corresponds to increasing of annealing time 
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Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 


Increasing of number of curve corresponds to increasing of annealing time 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and = y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and € = y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter & for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 


0.12 
0.08 
% 
Q 
ie) 
0.04 
0.00 
0.0 0.1 0.2 0.3 0.4 0.5 
a/L, &, & Y 


Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which have 
been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and &= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter & for a/L=1/2 and ¢ = y= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 
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Fig. 8. Normalized dependences of component uz of displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 


where y (x,y,z) is the approximation function. 
Dependences of optimal values of annealing time 
on parameters are presented on Figs. 6 and 7 for 
diffusion and ion types of doping, respectively. It 
should be noted, that it is necessary to anneal 
radiation defects after ion implantation. One could 
find spreading of concentration of distribution of 


dopant during this annealing. In the ideal case 
distribution of dopant achieves appropriate 
interfaces between materials of heterostructure 
during annealing of radiation defects. If dopant did 
not achieves any interfaces during annealing of 
radiation defects, it is practicably to additionally 
anneal the dopant. In this situation optimal value 
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of additional annealing time of implanted dopant 
is smaller, than annealing time of infused dopant. 


Further we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in 
doped areas of heterostructure. Under following 
condition «< 0 one can find compression of 
distribution of concentration of dopant near 
interface between materials of heterostructure. 
Contrary (at s>0) one can find spreading of 
distribution of concentration of dopant in this 
area. This changing of distribution of 
concentration of dopant could be at least partially 
compensated by using laser annealing [29]. This 
type of annealing gives us possibility to accelerate 
diffusion of dopant and other processes in 
annealed area due to inhomogenous distribution 
of temperature and Arrhenius law. Accounting 
relaxation of mismatch-induced stress in 
heterostructure could leads to changing of optimal 
values of annealing time. Mismatch-induced stress 
could be used to increase density of elements of 
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